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(57) Abstract 

A system and method for estimating Doppler shift compensation in a mobile communication system is described. A channel impulse 
response determination circuit determines a channel impulse response for the channel on which the signal is received. A reference vector is 
generated using a set of known symbols and the channel impulse response, and an error vector is then generated by comparing the reference 
vector with received samples of the signal. In addition, the noise variance of the incoming signal can be determined. Then, a Doppler 
correction factor is estimated having a real part representing amplitude correction and an imaginary part representing phase correction. 
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KSTTMATION OF DOPPLER SHIFT C OMPENSATION 
IN A MOBILE COMMUNICATION SYSTEM 

The present invention relates to estimating a Doppler shift 
compensation in a mobile communication system. 

In a mobile communication system, signals which are transmitted 
from mobile stations moving relative to a base station are 
subject to the well-known Doppler effect, which causes a 
frequency shift in the frequency received at the base station 
relative to that which was transmitted by the mobile station. 
This frequency shift is referred to herein as the Doppler shift. 
The Doppler shift is dependent upon the speed and direction of 
movement of the mobile station relative to the base station. 
Thus, the Doppler effect can provide an increase or a decrease 
in the frequency, depending on the direction of movement of the 
mobile station relative to the base station. The magnitude of 
the Doppler shift is dependent on the speed with which the mobile 
station is moving relative to the base station. 

Another area of frequency offsets arises when the transmitter 
unit and the receiver unit are wrongly synchronised. 

Existing mobile communication installations provide a form of 
Doppler compensation, in that the frequency detection circuitry 
within the base station which selects a particular signal on a 
particular channel can take into account a certain amount of 
Doppler shift in the signal. 

AU 664626 relates to a method and circuit arrangement for 
compensating for the Doppler shift in a radio signal propagating 
between a base station and a mobile station when the mobile 
station approaches and moves past the base station. As a mobile 
station approaches the base station with a decreasing propagating 
time delay, the propagating time delay is integrated at intervals 
to determine how the propagating time delay varies with time. 
This information is used to effect a change in the radio signal 
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frequency at a particular time to compensate for the sudden 
Doppler shift as the mobile station moves past the base station. 
In a TDMA system, the interrogation occurs at intervals equal to 
an integral number of time frames. In the GSM standard a time 
frame comprises eight consecutive time slots and a single 
transmission burst passes between a particular mobile station and 
base station in any one time frame. 

Thus, in this method of Doppler compensation decisions are made 
in response to past and incoming signals to improve the reception 
at a future time, i.e. it is a reactive system. It would be 
desirable to implement a system which can actively compensate for 
Doppler shifts in an incoming signal in real time. 

Moreover, due to fading a signal subject to the Doppler effect 
may have a variation in amplitude as well as frequency. It would 
be desirable to estimate the Doppler effect on amplitude as well 
as frequency. 

According to one aspect of the invention there is provided a 
system for generating a Doppler correction factor for 
compensating for Doppler shift in a signal transmitted between 
a mobile station and a base station in a mobile communication 
system, the system comprising: a channel impulse response 
determination circuit for determining a channel impulse response 
for the channel on which the signal is received; reference 
circuitry for generating a reference vector from a set of known 
symbols and the channel impulse response; means for generating 
an error vector from the reference vector and received samples 
of the signal; means for determining noise variance of the 
incoming signal; .and an estimator for generating a Doppler 
correction factor having a real part representing amplitude 
correction and an imaginary part representing phase correction, 
the estimator utilising the noise variance, the error vector, the 
samples of the received signal and the sampling times. 

According to the GSM standard, a TDMA transmission burst 
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comprises a training sequence which is normally used to determine 
the channel impulse response for the channel on which the signal 
is received. This is done by a convolution of the received 
training sequence with a stored version of the training sequence . 
Thus, the training sequence can constitute the set of known 
symbols for use in generating the reference vector. 
Alternatively, the set of known symbols can comprise estimated 
symbols from the received signal samples. 

In a TDMA system, the signal comprises a sequence of transmission 
bursts in respective time slots. The Doppler correction factor 
can be generated for each transmission burst and used to correct 
signal samples in that burst . 

In addition for correcting for the Doppler effect, the Doppler 
correction factor also corrects for frequency offsets arising 
from wrong frequency synchronisation between transmitter and 
receiver units . 

According to a further aspect of the present invention there is 
provided a method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: determining 
a channel impulse response for the channel on which the signal 
is received; generating a reference vector from a set of known 
symbols in the channel impulse response; generating an error 
vector from the reference vector and received samples of the 
signals; determining noise variance of the incoming signals; 
estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and applying the Doppler 
correction factor to the received signal samples thereby to 
compensate for Doppler shift . 

For a better understanding of the present invention and to show 
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how the same may be carried into effect refers* 
made by way of example to t-hm reference will now be 

Y or example to the accompanying drawings in which: 

Figure i is a diagram Qf a a 
communication system; mobile 

Figure 2 is a block diagram of circuitrv f n , < -, 
Doppler shift compensation; and ClrCUlt ^ f <~ elementing 

Figure 3 is a diagram illustrating the Estim fi block. 
Figure l illustrates a normal burst in a 

to th<= pc:m , , TDMA s y st em according 

signals IT ^ StaCi ° nS tr "~ lt bUrSCs »• -odul.J 

™ le ; n fr ;ri cy channels aiiocaced by * *»• 

oncroller. one frequency channel may support ur> fn • u„ 

: r r b c ssociaced with . im n 

call is allocated a time slot in which to send th- k 7 
Further details of a TDMA system according to the os M sta^d 

bits ,TS, are added at each end of ^ • tail 

training seguence ,«,, is a predetermined seguence of bits which 
is sent by the mobile station ,„s, and is kno™ at the base 

r ceiTd COntr ° ller <BSCK ^ Fi9Ure l ' " « that a 

received signal sample y k represents the ka bit of the ^ ^ J" 

this instance the start of the training seguence. The traininc 
sequence is utilise at- u training 

utilised at the base station controller to «v • 

the impulse response of the channel over which the burst TZT 

and at least in one embodiment of the invention to estimate a 

Doppler cessation factor. The actual information which is 

transmuted is located in the data bits (DATA) o, the burst 

As explained earlier, the environment through which a siona, 

passes from a mobile station to a base static 

case station can vary 
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considerably, depending, amongst other things, on the distance 
between the mobile station and the base station, and interference 
caused by buildings and other structures in the area. As a 
result, the signal strength and signal quality of the signal 
received at the base station varies widely. Moreover, for moving 
mobile stations, the signal received by the base station is 
subject to a Doppler shift which should be corrected. 

Figure 2 illustrates a circuit 1 suitable for implementing a 
Doppler compensation in a GSM system. It should be understood 
that the various blocks in Figure 2, although illustrated as 
separate interconnected entities, do not necessarily represent 
separate physical entities, but are intended to represent 
diagrammatically the various steps which are carried out. The 
blocks could be implemented as circuits or a suitably programmed 
microprocessor may effect each of the functions which is 
individually assigned to the blocks. 

An antenna 20 receives signals 11 from the mobile stations MS. 
The antenna 20 is connected to RF circuitry 22. This circuitry 
22 operates on the received burst to downshift the frequency to 
the baseband frequency and to sample the burst to provide from 
the analogue signal digital sampled values y t where i is the 
sample index. The RF circuitry 22 outputs the sampled values 
y 0 ...y n for each burst together with the sampling times t 0 ...t n , 
typically at the expected bit rate of the transmitted signal. 
The output of RF circuitry 22 is supplied along line 24 to a 
channel impulse response (C.I.R.) block 10, to a variance and 
error calculator 16 to filtering and equalisation circuitry 12, 
to an Estim /3 block 36 for estimating a Doppler correction factor 
and to transforming circuitry 4 0 to enable the application of the 
Doppler correction factor to the burst . 

A memory 32 holds the training sequence x which is the 
predetermined sequence of bits which is sent by the mobile 
station MS as a training sequence and received at the base 
station and converted into samples y^ = Y k -»y k425 . The reference 
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training sequence x is supplied to a reference generator 14 . „ 
to the channel impulse response (C.I.R., blocx xo Tn^eren 
generator 14 also receives the estimated channel impose " 
h fro, the channel i mpu i se response block ^ lmPUls * 

The c I.R. block 10 receives the burst, including the received 
traxn.ng sequence Ytrs and calculates an estimated channel ImpuLe 
response h for example by calculating the cross cox llZlll 
between the received training sequence y TOS and the Known t ^ 
sequence x according to equation 1. There are other type TZ 
estimation methods which could be used. 

h = xcorr (x,y TRS ) 

(equation 1) 

correlation ST""* *° e " eCti " 3 th * -oss 

correlation. the known training sequence x. which is stored in 

the training segu e„ ce has been modulated at the MS for 
transmission. according to the GSM standard. The cross 
correlation is done in a known manner to produce a chapel 
impulse response in the form of five tap values (h(j) . f 

AS is known, the estimated impulse response h is used to 
calculate the expected estimate of the data in the ITLZ 
burst, as though the data has heen subject to the same average 

The c I.R. block also generates timing advance information r 
which to deteraJne ^ ^ atio 

received burst is located. 

For each burst, the estimated channel impulse response h for that 
burst is calculated by the CIR block 10 and is supplied t o 

n ZTTT ng =irOUltry WhiCh all °" S "» — ' AT l<r 
in that burst to be recovered. As is known thl 

filtering/equalising circuit » receives the ch. "^ 
response h and timing information r for the received burst " 
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allow the signal to be demodulated, filtered and decoded, to 
recover the data in a known manner. 

The reference generator 14 produces a reference vector, z, which 
is calculated using the convolution of the impulse response and 
the known training sequence . The reference generator 14 performs 
the following calculation: 

z - h*x (equation 2) 

where x represents the known symbol values depending o the used 
modulation. 

z k represents the kth sample of the vector z and z denotes a 
diagonal matrix of the vector samples z k , that is, z=xh, where x 
is the diagonal matrix of the known training sequence symbols. 

The reference vector z is supplied from the reference generator 
to the variance and error calculator 16. As described above, the 
variance calculator also receives the samples y if including the 
received training sequence y^s- Ttle variance calculator 
calculates a variance var (a 2 ) according to the following 
equation: 

25 

var = E(|yi - z L | A 2)) (equation 4) 

i=4 



ref_length 

The term ref_ length is a constant representing the length of the 
reference vector, z. This is calculated by multiplying the 
number of samples (22) by the bit separation. 

In equation 4, the values of y A are the sampled values of the 
received training sequence for the burst (y k -*y k *2s) • 

The variance and error calculator also generates an error vector 
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z-y (i.e. error values z^yi for each sample). 

The error vector z-y and the variance a 2 are supplied to the 
Estim (3 block 36. The Estim fi block 36 also receives the signal 
samples and times y i# ti from the RF circuitry 22. The Estim /3 
block can be implemented as a processor executing a program or 
by other suitable means. It operates to generate a Doppler 
compensation factor 0 which is a complex number having real and 
imaginary parts. Thus, /S may be represented as: 

or+j^t 

where a represents amplitude variations due to fading and 0 is 
the Doppler angle. The Estim j8 block 3 6 also includes a memory 
for holding values required to generate the Doppler correction 
factor 0, including C 0I E(j8) and E(z-y). These terms are 
explained later with reference to Figure 3 . 

The Doppler correction factor /8 together with the signal samples 
Yt and times t t for each burst are supplied to the transforming 
circuitry 40 to implement Doppler correction according to the 
following equation: 

yn-compenaated ~ Yn^~^ n ) • 

Another implementation of the transforming circuitry 40 is to 
separate the amplitude and phase components into separate" 
correction steps. The phase can be corrected by multiplying the 
samples with w* n and the amplitude can be corrected ±n a 

separate step: 

yn-phase-compensated « y^'i'^W^ 

yn-compenaated = y-phaae-compenaated ( 1 " real ( j8) t„) (equation 5) 

According to the above, imag(/8)=0, and real (j8)=a. Thus, the 
received signal is either amplified or alternated during the 
burst depending whether real (0) is negative or positive. 
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This generates a set of compensated signal samples y comp which can 
be fed to the filtering/equalising circuitry 12 to generate 
refined values for the data of the burst • 

As an alternative to the above, in which the training sequence 
x is used to supply known symbols as a reference, bits estimated 
in an equalisation circuit 30 can be used. The equalisat ion 
circuit 3 0 , for example of Viterbi equaliser, receives the 
filtered, demodulated and equalised signal DATA (r) from the 
filtering and equalisation circuitry 12. The equalisation 
circuit 3 0 operates on a part of the data sequence data of the 
burst (that part having been derived from Estim. block in Figure 
1) to estimate output bits which were sent from the mobile 
station. This output is referred to herein as estim_bits . The 
equalisation circuit 30 operates to make decisions of the bits 
as in known mobile communication systems and thus it will not be 
described further herein. These estimated bits can be used in 
place of the training sequence to supply known values for 
calculating the reference vector z. The estim__bits are supplied 
to a reference circuit 34 with the channel impulse response h 
which generates an error vector z'-y, where z ' =estim_bits*h. 

The theory underlying the generation of the Doppler correction 
factor j3 will now be outlined. 

The problem of (pure) Doppler frequency estimation can be written 
as : 

Xh=Yr+w, where r is a rotation vector, 
Y is a diagonal sample matrix, and 
w is noise . 
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Xh corresponds to convolution between known symbols (diagonal 
-trxx x of training sequence symbols, and impulse response h 

So, in the Doppler frequency estimation, the DoppXer anole * 
to be estimated. T his problem is nonlinear and therefore a b t 
tr.cxy to be soived. particularly £ or indentation in rea 
trme. A linearisation step can be performed assuming e» % - 1+ H« t 
Due to fading the required Doppier frequency correction cannot 
be assumed to be constant in amplitude. 

If the term J*, has also a real component, the amplitude 
varrat.ons can roughly estimated. Thus , the Cerm P 
replaced by a more general assumption 1+ ,t„ where , is Comdex 

z i =y i (i+iet i )+w i , 

which can be reformatted in matrix notation as 
z-y=P/S+w 

(equation 6) 

where P=[y 0 t o y ltl . . .y n ,t n ,]T and fi ±g fco be estimated 

A solution to equation 6 can be found by using an LMMSE estimator 

where E(0) is an expected value for £ c r 
-variance matrixes between the two element, Zw'gX, Zl 

^■js".*^ ofT: n raft the noise is - hice *- 

y ariance ot (5 and after some manipulation 
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/3=E(i8)+p H ( Z -y-E(z-y) ) / (6 2 /C fi +p*p) , where 
C,«E(</8-E(jS) ) .(P-E(fi) ) H ) 



(equation 7) 
(equation 8) 



p H is the complex conjugate transpose matrix of p. 

Equation 7 is used in the Estim (3 block 3 6 to generate /? . The 
Estim /3 block is illustrated schematically in Figure 3. In 
practice, it can be implemented as software executed by a 
processor with a memory for. holding the required data. An 
Estimate_j3 block 50 operates on the data (Y^ti) received in each 
time slots to generate a value for (3 which is used by the 
transforming circuitry 40. An Acc (3 block 52 accumulates and 
averages (3 values over preceding time slots to generate E((3) . 
The Acc B function can be implemented by an adaptive filter. An 
Acc (z-y) block 54 generates an expected error E(z-y). To do 
this, it can use values for z-y from preceding time slots, or 
E(/3) from the Acc (3 block 52 to estimate a frequency error. A 
Gen C 0 block 56 produces a value for C 0 . A number of 
possibilities exist for the parameters of equation 7. 

If E((3) is set to zero, and 5 2 /C p is also set to zero, equation 
7 resolves to a normal linear regression. 

As mentioned above, it is intended that (3 will be generated for 
each time slot. The expected value E(/S) can be generated by 
calculating the average of estimated (3 values from previous time 
slots. Thus, these values can be retained as they are calculated 
for each time slot, and used to produce an average value held in 
the E((3) block in the Estim block 36. 

Another simplification is to assume that the real part of E(j3) 
is zero. This is a reasonable assumption because, in typical 
fading channel situations, no long time amplitude is expected. 
Thus, it is adequate to perform long time averages only in 
respect of error components of frequency. Also, the real part 
could be assumed to be zero if only phase variations are to be 
corrected. 
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It is possible to avoid having to calculate the expected error 
E(z-y) by derotating the received samples Yi by e l-«-»wt 
this purpose, fi can be an earlier estimated value of fi from °! 
previous time slot. In this case, equation 7 changes to 

MW + P«( 2 -y,/( jV C^, (eguat . on g) 

Thus, this requires a preprocessing step (see 50a in Figure 3) 
in which the sample values Yi . are derotated by Drior . 

forming the vector p " J 



The value for C (/>> can be generated by the GEN c, block using 
different approaches. It can be a constant value. I can be 
estimated using the information of previous time slots, using for 
example equation 8. It can be made a function of the expected 
value of fi and/or the noise variance 5*. it can be estimated 
using an adaptive filter, or any combination of the preceding. 

The technique is particularly suitable for short delay spread 
channels and pure frequency errors. 
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CLAIMS: 

1. A system for generating a Doppler correction factor for 
compensating for Doppler shift in a signal transmitted between 
a mobile station and a base station in a mobile communication 
system, the system comprising: 

a channel impulse response determination circuit for 
determining a channel impulse response for the channel on which 
the signal is received; 

reference circuitry for generating a reference vector from 
a set of known symbols and the channel impulse response; 

means for generating an error vector from the reference 
vector and received samples of the signal; 

means for determining noise variance of the incoming signal; 

and 

an estimator for generating a Doppler correction factor 
having a real part representing amplitude correction and an 
imaginary part representing phase correction, the estimator 
utilising the noise variance, the error vector, the samples of 
the received signal and the sampling times. 

2. A system according to claim 1, wherein the Doppler 
correction factor 0 is generated according to the following 
equation: 

0=E (£) +p H (z-y-E (z-y) ) / (SVC^p"?) 
wherein E(/3) is an expected value of /?, 

p H is a complex conjugate transpose matrix of p, where p is the 
vector (yi/ ti) T , 

z-y is the error vector, 

E(z-y) is the expected error in the received samples, 
6 2 is the noise variance, and 



JNSDOCID: <WO 0054431 A1_l_> 



BNS page 15 



WO 00/54431 



PCT/EP99/01618 



14 

C 0 is a constant. 

3. A system according to claim 2, wherein E(j8) and 5 a /C, are set 
to zero. 

4. A system according to claim 1, 2 or 3 , which comprises means 
for applying a Doppler correction to the received signal samples 
using the Doppler correction factor 0. 

5. A system according to any preceding claim wherein the signal 
comprises a transmission burst in a time slot in a TDMA mobile 
communication system, and wherein a new Doppler correction factor 
£ is generated for each time slot. 

6. A system according to claim 5, which comprises means for 
estimating E<0) as an average of generated values for fi from 
previous time slots. 

7. A system according to claim 5, wherein the real part of E(/8) 
'is set to be zero. 

8 - A system according to any preceding claim, which comprises 
means for derotating the received signal samples by e^^^^ 
where E(0) is an expected value for 0 the derotated received 
signal samples then being used to generate the Doppler correction 
factor fi in the present time slot according to the following 
equation: 

0=E(/J) +p H ( 2 -y) /(SVC^+p"?) . 

9. A method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: 

determining a channel impulse response for the channel on 
which the signal is received; 

generating a reference vector from a set of known symbols 
in the channel impulse response; 
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generating an error vector from the reference vector and 
received samples of the signals; 

determining noise variance of the incoming signals; 

estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and 

applying the Doppler correction factor to the received 
signal samples thereby to compensate for Doppler shift. 

10. A method according to claim 9, wherein the signal comprises 
a sequence of transmission bursts in a TDMA mobile communication 
system, and the Doppler correction factor is generated for each 
transmission burst. 

11. A method according to claim 9 or 10, wherein the Doppler 
correction factor j8 is generated according to the following 
equation: 

/3=E (0) +p H (z-y-E <z-y) ) / (fiVc^p"?) 
wherein E(jS) is an expected value of /?, 

p H is a complex conjugate transpose matrix of p, where p is the 
vector (y^tj*, 

z-y is the error vector, 

E(z-y) is the expected error in the received samples, 
5 2 is the noise variance, and 
is a constant. 
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